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We theoretically study slow collisions of NH3 molecules with He atoms, where we focus in partic- 
ular on the observation of scattering resonances. We calculate state-to-state integral and differential 
cross sections for collision energies ranging from 1CP 4 cm -1 to 130 cm" 1 , using fully converged 
quantum close-coupling calculations. To describe the interaction between the NH3 molecules and 
the He atoms, we present a four-dimensional potential energy surface, based on an accurate fit of 
4180 ab initio points. Prior to collision, we consider the ammonia molecules to be in their antisym- 
metric umbrella state with angular momentum j = 1 and projection k = 1, which is a suitable state 
for Stark deceleration. We find pronounced shape and Feshbach resonances, especially for inelastic 
collisions into the symmetric umbrella state with j = k = 1. We analyze the observed resonant 
structures in detail by looking at scattering wavefunctions, phase shifts, and lifetimes. Finally, we 
discuss the prospects for observing the predicted scattering resonances in future crossed molecular 
beam experiments with a Stark-decelerated NH3 beam. 



I. INTRODUCTION 

The tremendous experimental progress in performing 
scattering experiments has evolved to the point where it 
is nowadays possible to study collisions between particles 
in the laboratory over an energy range of about 25 orders 
of magnitude. The collisions of highest energy are pro- 
duced by modern charged-particle accelerators reaching 
the TeV range, while the collisions of lowest energy are 
studied in ultracold atomic quantum gases going all the 
way down to the nK regime. In both types of collision ex- 
periments scattering resonances play an important role. 
In high-energy collisions, a resonance in the cross section 
caused by the formation of an intermediate bound state is 
a direct way to detect previously unseen particles. In ul- 
tracold atomic scattering, the energy of an intermediate 
bound molecular state to be formed during the collision 
can sometimes be accurately tuned by applying an exter- 
nal magnetic field. As a result, the scattering length of 
low-energy s-wave collisions gets under full experimen- 
tal control, giving rise to a unique quantum many-body 
environment with a completely tunable interaction pa- 
rameter P, 

Because molecules arc typically harder to manipulate 
than atoms and charged particles, the observations of res- 
onances in molecular beam scattering have been limited 
to a few rare cases However, in recent years rapid 

progress has been made in performing high-precision cold 
molecular scattering experiments due to the application 
of the Stark deceleration technique to the study of molec- 
ular collisions @ . A Stark decelerator operates according 
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to the same principles as a linear charged-particle acceler- 
ator, where the dipolar or Stark force is used to decelerate 
neutral polar molecules with time-varying electric fields 
[i"o| . With the Stark decelerator it is possible to generate 
almost perfectly quantum-state selected molecular beams 
with a computer-controlled final velocity and a small lon- 
gitudinal velocity spread. By applying this technique to 
the scattering of the OH radical with rare gas atoms, 
such as Xe Q , Ar [ll[ , and He [l2| , the threshold behav- 
ior for inelastic scattering into the first excited rotational 
levels of OH could be accurately determined. Excellent 
agreement was found with cross sections obtained from 
close-coupling calculations using ab initio potential en- 
ergy surfaces (PESs) [{J [ll], [l3| . In the same way, also 
cold inelastic collisions of OH radicals with D2 molecules 
were studied experimentally 12 1. 

In this article, we study in detail cold collisions be- 
tween NH3 molecules and He atoms. The ammonia-He 
system is a van der Waals complex, and in general the 
(quasi-)bound states of such complexes are sensitive to 
the interaction potential. As a result, high-resolution 
spectroscopy on van der Waals complexes has been an 
important tool for increasing our understanding of inter- 
molecular forces [HI, EH- High-precision scattering ex- 
periments are a very promising additional tool for obtain- 
ing detailed information on potential energy surfaces. At 
higher scattering energies the short-range repulsive part 
of the interaction is mainly probed, while at very low col- 
lision energies the long-range part of the potential is dom- 
inant in determining the scattering behavior. Moreover, 
scattering resonances give important information on the 
energy of quasi-bound states that are sensitive to poten- 
tial wells at mid-range interparticle distances. This shows 
that large parts of the potential energy surfaces can be 
accurately probed by cold collision experiments. Recent 
scattering experiments have indeed been able to distin- 
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guish between PESs that were onl y of go od quality and 
PESs that were of excellent quality [ll|,[l3|. A very differ- 
ent experiment in which the NH3-HC interaction plays an 
important role, is the trapping of NH3 molecules inside 
He nanodroplets to perform high-resolution spectroscopy 

Rotational energy transfer by cold collisions is an im- 
portant process in various astrochemical environments, 
such as interstellar clouds and cold exoplanetary atmo- 
spheres. Since the first identification of NH3 molecules in 
the interstellar medium fl7j , ammonia has been detected 
in several gas-phase astrochemical spectra. The rate co- 
efficients of NH3-HC scattering arc an important ingredi- 
ent for a numerical modelling of astrochemical environ- 
ments. This is one of the reasons why NHs-He collisions 
have been studied experimentally flSl — f22j ] and theoreti- 
cally (23l - |28| by several groups. The most recent scatter- 
ing calculations have been performed with the potential 
energy surface of Hodges and Wheatley [29[. However, in 
order to get agreement with experimentally determined 
virial coefficients, this potential had to be scaled by a 
rather large factor |26|. The same potential has also 
been used to theoretically study low-energy NH3-He col- 
lisions, where strong scattering resonances were observed 
for various initial and final states of the NH3 molecule 
[2^ | . Unfortunately, the initial state that is most suit- 
able for Stark deceleration was not considered. This is 
namely the state \jk±) = |11— ), where j is the angular 
momentum of the ammonia molecule, k is the projec- 
tion on its threefold symmetry axis and +/— refers to 
its symmetric/antisymmetric umbrella inversion tunnel- 
ing state. For the energy level diagram of the ammonia 
molecule, see Fig. [TJ Moreover, in the study of Ref. [28| 
ammonia was treated as a rigid molecule, implying that 
the umbrella inversion motion of the NH3 molecule was 
not considered. 

In this article we study all possible elastic and in- 
elastic scattering processes at low collision energies, us- 
ing 1 11—) as an initial state of the para NH3 molecule. 
We show that particularly the inversion inelastic scatter- 
ing to the 1 11+) state gives rise to pronounced resonant 
structures that are promising to be observed experimen- 
tally in crossed beam experiments. We start by intro- 
ducing the theoretical framework for studying the atom- 
molecule collisions. After this, we present a new NH3- 
He potential using the most recent developments in elec- 
tronic structure calculations. We describe the numerical 
methods to fit the potential, after which we present the 
calculations of the integral and differential cross sections. 
In both cross sections we find rapid variations as a func- 
tion of energy, which are clear signs of resonant behavior. 
To determine the origin of these resonances we perform 
bound state calculations as well as reconstructions of the 
full scattering wavefunctions. The phase shifts and the 
lifetimes are also determined near resonance. Finally, we 
comment on the prospects of observing these scattering 
resonances in the NH3-HC system in the near future. 
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FIG. 1: Energy levels \jk±) of the ammonia molecule, where 
j is the angular momentum of the molecule, k is the projec- 
tion on its threefold symmetry axis and +/— refers to the 
symmetric/antisymmetric umbrella state. Throughout this 
article we use the 1 1 1 — ) state as the initial state, so that we 
only consider para ammonia. We do not take into account 
hyperfine interactions. The collision energy is defined relative 
to the initial state. 



II. THEORY 

To theoretically study the low-energy scattering of 
NH 3 molecules with He atoms we briefly introduce three 
coordinate frames that are used in the calculations [30j |. 
These coordinate frames are discussed in more detail in 
the Appendix. The first frame is an orthonormal, right- 
handed space-fixed ('sf') laboratory frame located at the 
center of mass Q of the dimer. The coordinate R is the 
length of the vector R that points from the center of 
mass X of the NH3 monomer to the He atom, while 8 sl 
is the zenith angle of the vector R and <fi sf is the az- 
imuth angle in the space-fixed frame. The second frame 
is an orthonormal, right-handed body-fixed ('bf') dimer 
frame, also centered at the center of mass of the dimer. 
As explained in the Appendix, this frame is obtained by 
a rotation that aligns its z axis with the vector R. The 
third frame is an orthonormal, right-handed monomer- 
fixed ('mf') frame centered at the center-of-mass of the 
NH3 molecule, whose z axis is aligned with the symme- 
try axis of the ammonia molecule. This monomer frame 
is obtained from the space-fixed frame by rotating over 
the three Euler angles ( si = (a si , /3 sf , 7 sf ). Here, a sf and 
/3 sf are the azimuth and zenith angles of the ammonia 
C3 symmetry axis in the space-fixed frame, while 7 sf de- 
scribes the rotation of the NH3 molecule about this axis. 
The Euler angles of the monomer frame can also be given 
with respect to the body-fixed frame, and are then de- 
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noted by £ = (a , /3 , 7 ). Finally, we introduce the 
umbrella or inversion angle p of the ammonia molecule, 
which is the angle between the z axis of the monomer 
frame and a vector pointing from the N atom to one of 
the H atoms, so that p — -k/2 corresponds to a planar 
NH3 geometry. 

The Hamiltonian of the NH3-He system can now be 
written as 
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(1) 



where throughout the article we set h = 1, -ff mo i is the 
Hamiltonian of the NH3 molecule, p is the reduced mass 
of the atom-molecule complex, j is the angular momen- 
tum operator of the NH3 monomer with respect to the 
body-fixed frame, J is the total angular momentum op- 
erator also with respect to the body-fixed frame, and Vint 
is the interaction energy. We consider the interaction po- 
tential to depend on four coordinates, which implies that 
we assume the N-H bond length to be fixed and NH3 
to keep its threefold symmetry. The Hamiltonian of the 
NH3 molecule includes the monomer's rotation, as well 
as the kinetic and potential energy of its umbrella motion 
[3l| . namely 
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where I xx (p), I yy (p) and I Z z{p) are the moments of in- 
ertia of the threefold symmetric ammonia molecule with 
respect to the monomer frame axes, see e.g. Ref. 3l| . 
while Ipp(p) = 3mH^o(cos 2 p + 77 sin 2 p) with toh the hy- 
drogen mass, ro = 1.9099ao the fixed N-H bond length 
[33 |. rj = mn/(3r?]H + ?tin) and 7tin the nitrogen mass. 
Moreover, we have that g(p) = I X xIyyIzzIpp, while the 
potential energy for the umbrella motion Vumb(p) leads 
to a double well potential that we model by 



Kmb(p) 



+ a p exp 



7T\ 2 
2 



(3) 



a p = 23 229 



with the parameters k p = 90 651 cm _1 rad~ 
cm -1 and b p = 3.1846 rad~ 2 . The resulting double well 
potential gives rise to umbrella vibration levels of which 
two levels have an energy below that of the barrier at the 
planar ammonia geometry |3lj . Moreover, each of these 
two vibration levels splits into a pair of energy levels due 
to tunneling. The parameters of the umbrella potential 
are chosen such that the experimental energy splitting 
between the two tunnel states in the vibrational ground 
state, as well as the experimental splittings [H[ between 
the ground state and the two tunnel states of the first 
vibrationally excited level are accurately reproduced. 

To treat the Schrodingcr equation in body-fixed coor- 
dinates, we expand the scattering wavefunction in the 



following coupled-channel basis 

* M (^) = ^£|n)x„(i?), 



(4) 



where the radial dependence of the wavefunction is given 
by Xn(-R), while the body- fixed angular basis set 



\j,k,K,j,Mjy 



(2j + l)(2J+l) 
32tt 3 
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is used to treat the angular part of the Hamiltonian. 
Here, I>££,(0 = e ima d^,(/3)e im ^ with d^,{p) the 
well-known Wigner d-functions, k is the projection quan- 
tum number of the monomer angular momentum with 
eigenvalue j on the monomer z axis, K is the projection 
quantum number of both the monomer angular momen- 
tum and the total angular momentum with eigenvalue J 
on the body-fixed dimer z axis, Mj is the projection of 
the total angular momentum on the space-fixed z axis, v 
is the umbrella vibration quantum number, and the su- 
perscript +/— refers to the even/odd umbrella tunneling 
function. 

As a result, our task is to solve the following second- 
order matrix differential equation 



d 2 X n'(R) 
dR 2 



= ]T(n'|W» Xn (i?), 



(6) 



where we introduced the operator W = 2p(E — H + K) 
with the kinetic energy operator K given by the second 
term on the right-hand side of Eq. ([T|) . We note that J 
and Mj are good quantum numbers, and that the oper- 
ator W is diagonal in J and independent of Mj. Fur- 
thermore, the monomer part of the Hamiltonian, H mo \, 
is also diagonal in the angular basis set. The complexity 
of the matrix (n'|VF|n) can be further reduced by consid- 
ering the symmetry properties of the NH3-He complex. 
Because the Hamiltonian commutes with permutations 
of the three hydrogen atoms in NH3 and the operator 
for inversion in space E* , it is useful to adapt the ba- 
sis states such that they transform as the irreducible 
representations of the corresponding molecular symme- 
try group Z?3/j (M) in the notation of Bunker and Jensen 
34j . The adapted basis states of different symmetry can- 
not be mixed by the Hamiltonian. The precise procedure 
for this adaptation is described in the Appendix. More- 
over, in Ref. [l4| several useful relations can be found for 
determining the matrix elements of the W operator in 
the angular basis. 

In order to fully solve Eq. © , the wavefunctions must 
satisfy the appropriate scattering boundary conditions 
[35j . These boundary conditions are directly formulated 
in a space-fixed frame. The exact solution of the space- 
fixed Schrodingcr equation at larger separations R, i.e., 
when the interaction energy has approached zero, is a lin- 
ear combination of the proper spherical Bcsscl functions. 
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These Bessel functions are labelled by the space-fixed 
end-over-end rotational quantum number L, which has 
become a good quantum number at such large separa- 
tions. Therefore, the matching of the propagated wave- 
function from Eq. (J6j) to spherical Bessel functions can be 
performed at distances where the centrifugal energy, set 
by L and decaying as R~ 2 , is still large, as long as the in- 
teraction energy, decaying in our case as i? -6 , has become 
negligibly small. To perform the matching, it is necessary 
to transform between the body-fixed and the space-fixed 
basis sets. The latter basis set is for the present case 
given by 



\j,k,L, J,Mj,v±) 



(2j + l)(2i + l) 
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Mr. 



with Cm l (8 s{ , 4> si ) the Racah- normalized spherical har- 
monics and Ml the projection of the end-over-end angu- 
lar momentum on the space-fixed z axis, uij the projec- 
tion of the monomer angular momentum on the space- 
fixed z axis, and (jirni; J2m2|j3m3) a Clcbsch-Gordan 
coefficient. The transformation between the body-fixed 
and the space-fixed basis then becomes [l4[ 



\j,k,L,J,Mj,v±) (8) 

/or. _l i \ V2 

= J2\j^,K,J,M J ,v ± ) ^— T J (jK;L0\JK). 



To end our discussion of the matching procedure, we 
mention the various possible open and closed channels 
following from the Hamiltonian of Eq. @. The vibration- 
tunneling states of the umbrella motion are determined 
by calculating the cigenstatcs of the Hamiltonian formed 
by the last two terms of Eq. ([2]). Only the lowest 
four eigenstates, labeled (j)t(p) with vibrational quan- 
tum numbers v = and v = 1, are kept. With these 
four states as a basis for the umbrella motion, we turn 
to the first term of the Hamiltonian of Eq. $Z§. As a 
result, the rotational constants \/2I\\(p) become 4x4 
matrices, but the + and — states are not mixed. The 
eigenstates that result from diagonalization of the full 
monomer Hamiltonian in a basis containing all rotational 
states with j < 6 and the four umbrella states, deter- 
mine the open and closed channels. We label the open 
channels by \jk±); the vibrational quantum number v is 
omitted from this notation because all vibrational states 
with v > are closed for the energy range in which we 
are interested. We consider the ammonia molecules to be 
prepared in the |11— ) state, so that the lower lying |11+) 
state is open for all collision energies. Increasing the colli- 
sion energy beyond the energy of excited monomer states 
opens up the corresponding channels, and inelastic scat- 
tering into these states occurs if it is allowed by symme- 
try. The matching procedure of the wavefunctions to the 
boundary conditions for scattering at large distance R 



ultimately leads to an expression for the scattering ma- 
trix. This S matrix is subsequently directly related to 
the differential and integral cross sections for the elastic 
and inelastic channels, which can be compared with the 
outcome of collision experiments [36| . 



III. THE NH3-HE POTENTIAL 

Before we can apply the above described formalism 
to solve the scattering problem, we need to determine 
the NH3-He interaction potential. To this end ab initio 
calculations were performed with molpro [13 , using the 
supcrmolcculc approach with the counterpoise procedure 
of Boys and Bernardi [38[ . We considered the interaction 
energy to be dependent on four coordinates, namely R, 
/3 bf , 7 bf and p. 

The grid for the ab initio calculations consisted of 4180 
points. For R, in total 19 points were used. In the short 
and intermediate range, i.e., for 4eto < R < lOeto, we 
used an equidistant grid of 13 points with a separation 
of 0.5ao, while in the long range, that is for R > 10ao, 
we used an approximately logarithmic grid consisting of 
the points 12ao, 14.4ao, 17.3ao, 20.8ao, 25ao, and 30ao. 
For j3 , we used a Gauss-Legendre grid consisting of 11 
points for < /3 bf < tt, while for 7 bf , we used an equidis- 
tant Gauss-Chebyshev grid consisting of the points 7r/24, 
37r/24, 57I-/24, and 77r/24. Finally, for the grid in p we 
used an equidistant grid of 5 points, where the middle 
point was given by the value p% = 0.62267T, while the dis- 
tance between the points was given by Ap = (2p3 — 7r)/5. 

The calculations in the long range were performed 
with the coupled-cluster method taking into account sin- 
gle and double excitations and a perturbative treatment 
of triple excitations [CCSD(T)], using the augmented 
correlation-consistent polarized valence quadruplc-zeta 
(AVQZ) basis set. For the short range we used the ex- 
plicitly correlated CCSD(T)-F12 method [39| to account 
more efficiently for the strong effect of electron corre- 
lations in this regime. The CCSD(T)-F12 method was 
found to yield accurate results with the smaller triplc- 
zeta basis set (AVTZ), as illustrated by Table |TJ In 
this table the interaction energies are shown for differ- 
ent NH3-He geometries in the short, intermediate and 
long range with both the CCSD(T) and the CCSD(T)- 
F12 method using different basis sets. Also the effect 
of using midbond functions is included in this table, 
where the midbond orbitals were located along the vector 
connecting the center of mass of the ammonia molecule 
with the helium atom at a distance of (ro+R)/2 from the 
ammonia center of mass. Particularly in the short and 
intermediate range these midbond functions improve the 
interaction energies, so that they were used in the cal- 
culations with F12. From Table [H we see that in the 
short range the CCSD(T)-F12 method with an AVTZ 
basis set including midbond functions performs better 
than the CCSD(T) method with an AV5Z basis set, al- 
though the latter calculation is much more expensive due 



5 



TABLE I: Comparison of the interaction energy between the 
CCSD(T) method (abbreviated as CC) and the CCSD(T)- 
F12 method (abbreviated as F12) for different basis sets and 
different geometries as calculated with molpro [37j |. We used 
the augmented triple zeta (AVTZ), quadruple zeta (AVQZ) 
and quintuple zeta (AV5Z) basis sets. We also studied the 
effect of midbond functions which are indicated in the 
Table by the + sign, when they are added to the basis set. 
For the short-range geometry, indicated by x s in the Table, 
we used R = 4.5a , /3 sf = 0, 7 sf = and p = 14tt/24. For 
the mid-range geometry, indicated by x m , we used R = 7oo, 
/3 sf = 7r/2, 7 sf = 7r/6 and p = 157r/24. For the long-range 
geometry, indicated by xi, we used R = 15ao, /3 sf = 7T, 7 sf = 
7r/3 and p — 167r/24. The interaction energies are given in 
cm -1 . 

AVTZ AVTZ+ AVQZ AVQZ+ AV5Z 

CC (x s ) f446.2 1414.2 1407.4 1397.3 1394.7 
F12 (xs) 1393.5 1386.0 1389.8 1386.8 1386.1 
CC (x m ) -21.716 -23.529 -22.733 -23.521 -23.179 
F12 (x m ) -23.329 -23.679 -23.237 -23.510 -23.381 
CC (xi) -0.2506 -0.2583 -0.2521 -0.2538 -0.2526 



F12 (xi) 



-0.2622 -0.2637 -0.2548 -0.2554 -0.2524 



to the large basis set. 

The main reason why we did not use the F12 method in 
the long range is that we found that the implementation 
of this method in molpro (37} gives rise to an incorrect 
1/R behavior in the very long range, rather than the cor- 
rect 1/R 6 behavior for the system under consideration. 
This behavior is caused by the fitting of the electron den- 
sity distributions, which unfortunately does not result in 
exactly charge neutral monomers. Although the artifical 
residual charges can be reduced by introducing a larger 
electron density fitting basis set, the 1/R behavior will 
eventually always dominate the correct 1/R 6 behavior. 
Hence, we decided to use the CCSD(T) method without 
F12 for the long range. We used the AVQZ basis set, 
and we may conclude from Table U that this basis set 
indeed gives rise to accurate interaction energies in the 
long range. In order to switch smoothly between the re- 
sults of the two methods, we used the switching function 
s(R) 



s(R) 




if R<a 
if R > b 

\ sin (3 — sin 2 ^ ) otherwise 



(9) 



with x = (2R— b— a)/(b— a), a = 10ao and b = 13arj. The 
function is chosen such that the first three derivatives of 
s at R = a and R = b are zero. We thus calculated the in- 
teraction energies for the angular geometries at distance 
R = 12ao with both methods, where the calculated value 
with F12 was given a weight of 1 — s(12ao), while the 
value without F12 was given a weight s(12ao). 

To obtain an analytic representation of the interaction 
potential between the NH 3 molecules and the He atoms 
we first perform an expansion in tesseral spherical har- 
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FIG. 2: Contour plot of the NUVHe interaction potential 
as a function of R and /3 bf for 7 bf = and the equilibrium 
umbrella angle p c = 112.15°. The energies of the contours are 
given in cm" 1 . 



monies, namely 



^ nt (i?,/3 M ,7 M 



/ . i n 



( (i?,p)S ;m (/3 bt , 7 bt ) 



(10) 

where due to the symmetry of the dimer only terms with 
m = 0,3,6,... are present. Because we have 11 grid 
points in j3 hi , the summation over I is from < I < 10. 
The summation over m is from to the largest multiple 
of 3 that is smaller than or equal to the corresponding / 
value. On all grid points Ri and pj we determine the co- 
efficients of the angular expansion vi m (Ri, pj) by means 
of a quadrature on the ab initio grid with the ap prop riate 
Gauss-Lcgcndrc and Gauss-Chcbyshev weig hts hj. For 
the resulting expansion coefficients U; m , we distinguish 
between the short-range and the long-range behavior, so 
that vim (R,p)= vfm (R, p) + vf m (R,p). Both in the short 
range and the long range the dependence of the coeffi- 
cients vi m on p is represented by a polynomial expansion 
in (p — 7r/2) p , where p ranges from to 9. If I + m is 
even, then the polynomial expansion only contains even 
powers in p, while if I + m is odd only odd powers are 
present. In the long range, we expanded the potential in 
powers of R~ n , resulting in 



vf m {R,p) 



cimpnMaR) (p - |) P R- n . (11) 

n,p 



where the inverse powers of R that are involved depend 
on I. It can be shown [4l| that for I = 0, 2 the expansion 
starts with ni = 6, while for I = 1, 3 it starts with rij = 7. 
For I > 4, it starts with = I + 4. We used the analytic 
long range expansion of Eq. ([HI) only for I < 5, and for 
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each I we took the leading term R rii and the next-to- 



lcading term R n - 
damping function 



into account. The Tang-Toennies 



/»(») = !" E| 



(12) 



,i=0 



was included to avoid the singular behavior of the 
long-range terms in the short range j42j . For the 
value of a in Eq. we used the isotropic expo- 

nent in the short-range, or, to be more precise, a = 
ln[voo(Ri,P3)/v 00 (R 2 ,p 3 )}/AR = 2.088ag 1 with R x and 
i?2 the first two points of the R grid and AR = R 2 — Ri . 
The expansion coefficients c; TOpn were obtained from 
vim{Ri, Pj) by performing a weighted least-squares fit us- 
ing the last three points of the R grid and all points of 
the p grid. The three R points were weighted for each 
I by with m the leading power of the long range 
decay for the considered I. 

To describe the short and intermediate range of the po- 
tential the same expansion was employed in /3 bf ,7 bf and 
p as for the long range, but the behavior in R was treated 
differently. The following procedure was used. First, for 
all the grid geometries, the corresponding value of the 
analytic long range potential was subtracted from the ab 
initio values. Next, after performing the expansion in 
tesseral harmonics and powers of (p — 7r/2), the behav- 
ior of the resulting coefficients vi mp {R) was interpolated 
with a reproducing kernel Hilbert space (RKHS) method 
[43j . The smoothness parameter of the RKHS interpola- 
tion was set to 2, while the RKHS parameter to, which 
determines the power with which the interpolation func- 
tion decays, was chosen to depend on I. For I < 5, the 
parameter was set to m = ni + 1. Then the RKHS func- 
tion decays as R~ ni ~ 2 , which is faster for each I than the 
leading term in the analytic fit of Eq. (TTTj) . However, for 
I > 5, no analytic long range fit was done, and we used 
to = n, — 1. As a result, the corresponding RKHS func- 
tions decayed for each I as R~ ni with m — I + 4, which 
is the correct leading long range behavior for I > 5 [4l| . 

We have compared the fitted potential with the ab ini- 
tio values on the full grid to test the accuracy of the fit 
in the angles /3 bf ,7 bf and p. The quality of the fit in R 
cannot be tested in this way, because the RKHS proce- 
dure goes by construction precisely through the points 
to be fitted. We calculated the RMS (root mean square) 
error for each grid distance Ri and divided by the mean 
ab initio interaction energy at that distance, giving for 
the relative RMS error £(i?i) that 
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V^ bl (i?i, /3 bf , 7 bf , pi), and the summations are over all 
n = 220 angular grid points. For our potential fit we 
found that the relative error £(Ri) is less than 0.05 % for 




FIG. 3: Coefficients vi m (R,Pe) of the NH3-He interaction en- 
ergy as a function of center-of-mass distance R, evaluated at 
the equilibrium umbrella angle p e = 112.15°. The isotropic 
voo(R, Pe) coefficient is largest. The W10, t>20, ^30, V33, U40 and 
W43 coeffients are shown as well. 



all Ri, so the fits in the angular coordinates are excellent. 
To also test the fit in R, we calculated ab initio interac- 
tion energies for an additional 495 points, that were cho- 
sen to lie about halfway between the grid points used for 
the fit. The relative RMS error of the values calculated 
from the fit compared to the new ab initio values de- 
pended quite strongly on R, where the largest error was 
found to occur in the short range. Namely, for the test 
points at R = 4.3ao we found with the use of 45 different 
angular points a relative RMS error of 3.5%, while for all 
other R values we obtained a relative error of about 0.5% 
or less. An important reason for this behavior is that we 
use a RKHS fit for the short range, which behaves as 
a power law, while the true behavior of the potential is 
exponential. The fitting procedure could thus have been 
further improved using an exponential form. However, 
we note that already the present fitting error is rather 
small. Moreover, in the present paper we use the po- 
tential to describe cold scattering with collision energies 
of maximally 130 cm -1 , so that the extreme short-range 
behavior of the potential is not being probed. 

In Fig. [5J we show a contour plot of the NH3-He inter- 
action potential for 7 bf = and the equilibrium umbrella 
angle p e = 112.15°. For this value of 7 and p, the min- 
imum of the potential is given by V nlln = —35.08 cm -1 
for R = 6.095 a and /3 b£ = 89.0°. This may be com- 
pared to the potential of Hodges et al. (29|, where the 
minimum of the potential for 7 bf = and p e = 112.15° 
is given by V mnl = —33.46 cm" 1 for R = 6.133 ao and 
(3 hl = 88.75°. Although this difference in the well depth 
is not very large, we have found that the consequences for 
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low-energy scattering can still be quite substantial, as we 
will discuss in Section [V] Finally, we have for the leading 
isotropic coefficient, defined as Cq = ~v 0Q (R, p 3 )R 6 for 
large R, that in atomic units Cq = 39.6 Eh ■ The 
relative importance of the various vi m {R, p e ) expansion 
coefficients is shown in Fig. Ep The potential is available 
in FORTRAN 77 on epaps |44|. 



IV. COMPUTATIONAL ASPECTS 

Having discussed the formalism and the potential, we 
now turn to the numerical procedures that we used in 
order to obtain converged cross sections that can be com- 
pared with future cold-collision experiments. To numer- 
ically determine the four lowest lying vibration-inversion 
levels 4>v{p) of the Hamiltonian of Eq. for j = 0, we 
used the discrete variable representation based on sinc- 
functions (sinc-DVR) [ij|. The resulting eigenfunctions 
were used to determine the matrix elements (<fi'\l/2I\\\(f>) 
with (j) = <f>y{p) by numerical integration. For the prop- 
agation in solving Eq. ((HJ), the renormalized Numerov 
algorithm was used, starting at 4an and ending at 50ao, 
using an equidistant grid with 273 points. The renormal- 
ized Numerov method also allows for a complete recon- 
struction of the scattering wavefunctions. 

The angular basis set contained all monomer states 
with j < 6, where we checked that the inclusion of more 
monomer levels resulted only in deviations of maximally 
1 % for the calculated cross sections. The maximal value 
for the total angular momentum J that we used depended 
on the collision energy. For collision energies E < 10 
cm -1 , we included all angular basis states with J < 10, 
while for 10 < E < 50 cm , we included all basis states 
with J < 20, and for 50 < E < 130 cm -1 , we included 
all states with J < 30. The convergence of the cross 
sections with respect to the total angular momentum J 
is slowest for the elastic cross section. The inelastic cross 
sections are converged at considerably lower values of J 
than reported here. 

In order to check our results and gain additional in- 
sight, we also implemented a commonly applied model 
to treat the ammonia umbrella motion in scattering cal- 
culations [46|, [43]. In this model, no vibrationally ex- 
cited umbrella states are taken into account and the 
ground-state umbrella tunneling states are approximated 
as an even and odd combination of the two rigid equilib- 
rium structures. These two states are thus written as 
l±> = U(P - Pe) ±f(7T-p + p e )]/2 1 /2, W hcre f(x) is 
a function localized around x — 0. More precisely, the 
two-state model amounts to approximating the matrix el- 
ements of the potential by (±|v; m (i?, p)|±) = Vi m (R,p e ) 
for I + m even, and (±\vi m (R, /o)|=f) = Vi m (R,p e ) for 
/ + m odd. For the rotational constants we use the ex- 
perimentally determined values A xx = A yy = 9.9402 
cm' 1 , and A zz = 6.3044 cm -1 in the model. Further- 
more, we include the experimental ground state splitting 
of 0.79 cm -1 [33|] between the two tunneling states in the 
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FIG. 4: Integral cross sections for NH3-He scattering as a 
function of collision energy. The initial state of the NH3 
molecule is the |11— ) state. At low collision energies only 
elastic scattering (upper red curve) and inelastic scattering 
into the lower lying |11+) state (lower blue curve) can occur. 



scattering calculations. This simple model has been im- 
plemented in the scattering program MOLSCAT [4^]. We 
have used MOLSCAT to double-check the results that we 
obtained from our own scattering program. The model 
was previously found to result in good agreement with 
more elaborate treatments of the umbrella motion for 
scattering at higher collision energies (49| . In this article 
we also want to test the accuracy of the model for cold 
collisions, and in particular for the calculation of scatter- 
ing resonances. 



V. RESULTS 

In Fig. @J we show the integral cross sections for the 
scattering of NH3 molecules with He atoms for collision 
energies ranging from 10 -4 cm^ 1 to 20 cm' 1 . Initially, 
we only have elastic scattering and inelastic scattering 
into the |11+) state, which lies 0.79 cm -1 lower in energy. 
Figure S] was made using the previously described elabo- 
rate treatment of the umbrella motion, however, with the 
use of the model treatment almost exactly the same re- 
sults were obtained. We observe in the first place that, in 
agreement with the Wigner threshold laws [50|, the elas- 
tic cross section becomes constant for very small collision 
energies, while the inelastic cross section decreases with 
E as 1/yE. Going more into the details of the figure, we 
observe two shape resonances closely together in the elas- 
tic channel at collision energies of 1.86 and 2.22 cm -1 . 
In bound state calculations with the NH 3 -He complex 
enclosed in a box of variable size we found continuum 
levels with nearly the same energies that are practically 
independent of the box size, so we may conclude that 
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FIG. 5: Elastic (upper red curve) and inversion inelastic 
(lower blue curve) integral cross sections as a function of col- 
lision energy for NH3-He scattering. The inelastic cross sec- 
tion is scaled with a factor of 150, so that the actual inelastic 
cross section is 150 times smaller than shown in the figure. 
At higher collision energies the |22±) and the |21±) channels 
open. As a result, Feshbach resonances are observed, which 
are most pronounced for the inelastic scattering into the |11+) 
state. 



these peaks in the scattering cross section indeed corre- 
spond to shape or orbiting resonances caused by quasi- 
bound states. Such quasi-bound states may occur either 
in the incoming or in the outgoing scattering channel; for 
the specific case of elastic scattering these are the same. 
Looking at the dominant contributions to the cross sec- 
tion, the first peak was found to be mainly caused by 
quasi-bound states with total angular momenta J = 4 
and J = 6, while the second peak was mainly caused by 
a quasi-bound state with total angular momentum J = 5. 
In both cases they corresponded to an end-over-end an- 
gular momentum of L = 5. Looking in the same energy 
range at the inelastic scattering into the |11+) state, we 
observe not only two similar peaks at the same collision 
energies, but also two additional peaks at 1.08 and 1.44 
cm -1 . These two additional shape resonances can be 
readily understood by noting that for inelastic scattering 
the resonant quasi-bound state can occur either in the 
incoming channel or in the outgoing channel, where the 
latter channel is about 0.8 enr 1 lower in energy. This 
is indeed precisely the energy with which the two addi- 
tional peaks in the inelastic channel are shifted to the left 
in Fig. [U 

For completeness, we note that we studied this colli- 
sion energ y ra nge also with the potential of Hodges and 
Wheatley [29[ . Although the difference in the well depth 
between the two potentials at the equilibrium umbrella 
angle was only about 5%, we still found large differences 
in the resonant structures at very low energies. For ex- 
ample, using the Hodges and Wheatley potential [29| . we 
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FIG. 6: (a) Contributions of various channels and partial 
waves to the scattering wavefunction at a collision energy of 
37.28 cm -1 , where a Feshbach resonance for inelastic scat- 
tering into the |11+) state occurs. We used a total angular 
momentum of J = 3 and considered E" symmetry. As a 
result, the initial state 1 1 1 — ) asymptotically corresponds to 
L = 3, while there are eight open partial waves in outgoing 
channels, (b) Contributions of various channels to the square 
of the wavefunction. In the inner region, where the collision 
takes place, a large amplitude in the asymptotically closed 
1 21 — ) channel is observed. 



observed two very strong shape resonances at collision 
energies of 0.03 and 0.45 cm -1 induced by quasi-bound 
states with total angular momenta J = 3 and J = 4 
and end-over-end angular momentum L = 4. However, 
because our own potential is deeper, we find that these 
quasi-bound states have become true bound states with 
energies below the scattering continuum, so that they 
cannot cause shape resonances anymore. As a result, 
the first shape resonances we find with our potential are 
induced by quasi-bound states with total angular mo- 
menta J = 4, 5 and 6 and L = 5, as shown in Fig. 2) 
This point also clearly shows that scattering resonances 
at low energy can be very sensitive to the precise shape 
of the potential energy surface, which means that accu- 
rate scattering experiments can be used to probe very 
precisely our knowledge of intermolecular interactions. 

In Fig. [5j we again show the integral cross sections for 
elastic scattering and inelastic scattering into the |11+) 
state, but now considering collision energies from 10 to 50 
cm^ 1 . Note that the inelastic cross is actually 150 times 
smaller than shown in the figure. As can be seen from 
Fig.[U at a collision energy of 28.33 cm -1 it becomes en- 
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FIG. 7: Top panel: Phase shift sum as a function of collision 
energy for various total angular momenta J. Both E' and 
E" symmetries are considered, so that we have two curves 
for each J. Bottom panel: the corresponding lifetimes as a 
function of the collision energy. The lifetimes are obtained 
from the derivative of the phase shifts with respect to the 
collision energy (in cm -1 ); they are given in units of 5.3088 
ps. 



ergetically possible to excite the ammonia molecule from 
its 1 1 1 — ) state to its |22+) state, and a new scattering 
channel opens. At 29.12, 39.33 and 40.12 cm -1 , the 
[22 — ), |21+) and |21— ) channels open, respectively, as 
also indicated in Fig. [5] The opening of the new channels 
is seen to have a profound effect on the inelastic cross sec- 
tions to the 1 11+) state. Namely, before these new chan- 
nels open a bunch of Feshbach resonances is observed. 
These resonant structures are called Feshbach resonances 
because they are caused by a molecular level that is dif- 
ferent from the incoming and the outgoing channel. In 
Fig. [5] we see that especially the Feshbach resonances in- 
duced by the |21±) levels at collision energies around 
40 cm -1 are strong, giving rise to almost a factor of 3 
increase compared to the background inelastic cross sec- 
tion. These resonances seem to be particularly suited to 
observe in a collision experiment. We come back to this 
point more elaborately in Section [VII 

To understand the Feshbach resonances in more de- 
tail, we have studied the scattering wavefunctions. In 
Fig. |5] contributions to the scattering wavef unction are 
shown at a collision energy of 37.28 cm -1 . At this col- 
lision energy, there is a Feshbach resonance for inelastic 
scattering into the |11+) state, somewhat below the en- 
ergy at which the |21— ) channel opens. In Fig. HJa), we 
show contributions of different open and closed channels 
to the scattering wavefunction. For this particular fig- 



ure, we considered a total angular momentum of J = 3 
and symmetry E" (see the Appendix). This means that 
for the incoming channel, i.e., the |11— ) state, asymp- 
totically only the partial wave with L = 3 contributed. 
For the four open outgoing channels, namely 1 1 1=1=) and 
|22±), in total eight open partial waves are possible for 
the considered J and E" symmetry. In the inner re- 
gion also contributions corresponding to asymptotically 
closed channels can gain amplitude, when they are cou- 
pled to the considered incoming state and outgoing state 
by the interaction potential. In Fig. EJJb), we show for 
each channel the resulting contributions to the square of 
the wavefunction. From Figs. [Ha) and (b) we clearly 
see that in particular the closed |21— ) channel has a very 
strong amplitude in the collision region, which shows that 
this state is responsible for the strong Feshbach resonance 
observed at this collision energy. 

A different way to study the Feshbach resonances is 
by looking at phase shifts in the scattering wavefunction. 
These phase shifts can be obtained from the eigenvalues 
of the scattering matrix [36|, [5l|. In Fig. [71 we show in 
the top panel the sums of the phase shifts in all open 
channels for various total angular momenta J. Since we 
consider both symmetries E' and E" , we have two curves 
for each J. From scattering theory it follows that when 
a resonance occurs, the phase shift sum rapidly increases 
by 7r [5l| as a function of energy. In the top panel of 
Fig. [Jj we indeed see this happening many times at the 
collision energies where resonances are found in the elas- 
tic and inelastic cross sections. The derivatives of the 
phase shifts with respect to the energy give the lifetime 
of the collision complex [36| . These lifetimes are shown 
in the lower panel of Fig. [7] This figure shows that at 
the collision energies where resonances occur, we indeed 
have long-lived intermediate quasi-bound states. 

In Fig. we show the integral cross sections for in- 
clastic scattering into the |21±) and the |22±) states, 
for a collision energy ranging from the energies at which 
these channels open, up to 50 cm -1 . Note that the in- 
elastic cross sections to the |22+) and |21+) states are 
scaled with a factor of 2. Immediately after each channel 
opens, we see strong resonant features, which are shape 
resonances, caused by quasi-bound states in the outgo- 
ing channel. In the |22±) channels we also find Feshbach 
resonances due to quasi-bound states of |21±) character. 

In the energy range from 10 to 50 cm -1 , we also stud- 
ied the scattering cross sections using the previously de- 
scribed model treatment of the NH3 umbrella motion. 
We found that the model calculations have the tendency 
to somewhat overestimate the strength of certain reso- 
nance peaks compared to the elaborate treatment of the 
umbrella motion. Studying this effect in more detail, we 
found that the differences are mainly due to the approx- 
imation of the nonzero potential matrix elements for the 
two tunneling states as vi m (R, p e ), rather than due to the 
neglect of the higher lying (j)f(p) states. Namely, by cal- 
culating the cross sections with the elaborate treatment 
and taking only the lowest two umbrella functions cf>^ (p) 
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FIG. 8: Inelastic integral cross sections for scattering into 
the |22±) channels and the |21±) channels as a function of 
collision energy. The inelastic cross sections for the |22+) 
and the |21+) channels are scaled with a factor of 2. After 
the various channels open, shape resonances are observed. 
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FIG. 9: Elastic and inelastic integral cross sections for scat- 
tering into the |11+), |22±), |21±) states as a function of 
collision energy. The inelastic cross sections for the |22— ) 
and 1 21—) channels are scaled with a factor of 20, and for 
the |11+), 1 22+) and |21+) channels with a factor of 40. At 
higher collision energies various |3fc±) and |4fc±) channels for 
the para ammonia molecules open. As a result, small Fesh- 
bach resonances are observed. 



into account we obtained cross sections that were nearly 
equal to the elaborate treatment with four umbrella func- 
tions, while they gave rise to the same differences with the 
model treatment. However, we note that in general the 
model treatment performed very satisfactory in describ- 
ing the resonance structures. All resonant peaks found 
with the elaborate treatment were also found with the 
model treatment, and typically the strength of the scat- 
tering resonances differed by less than 10%. Because the 
precise strength and location of the resonances are very 
sensitive to the the potential, we conclude that the use of 
the model treatment is useful in studying scattering reso- 
nances, especially in cases when the elaborate treatment 
is computationally too expensive. 

Looking again at Fig. [51 we note that there are signif- 
icant differences in the magnitudes of the inelastic cross 
sections for the various collision channels. For exam- 
ple, the transition to the |22— ) state is seen to be much 
stronger than the transition to the |22+) state, and the 
same holds for the transition to the |21— ) state compared 
to the |21+) state. The relative magnitude of the inte- 
gral cross sections for the elastic and inelastic scattering 
channels can even be more clearly observed in Fig. [SJ In 
this figure, we show the integral cross sections for scatter- 
ing into the |11±), |21±) and |22±) states, for collision 
energies ranging from 80 cm" 1 to 130 cm -1 . Notice the 
scaling of the inelastic cross sections indicated in the fig- 
ure. To explain the relative strengths of the transitions 
shown, we note that the scattering from the |11— ) chan- 
nel into different \jk±) channels is caused by different 
anisotropic terms in the interaction potential with coef- 
ficients vi m (R,p), cf. Eq. (JTDJ). For example, in order to 



change the umbrella state of the ammonia molecule (i.e., 
going from the odd — state to the even + state) we need 
terms in the potential for which I + m is odd, so that 
also the corresponding coefficient vi m (R, p) is odd in p. 
The various potential energy coefficients are plotted as 
a function of R at the equilibrium umbrella angle p e in 
Fig. El 

From this figure, we see that the isotropic coefficient 
vqo(R, p e ) is by far the largest coefficient of all. This coef- 
ficient can only couple the initial 1 1 1 — ) state to itself (see 
for example Refs. [H|,[52|]), causing a large elastic scatter- 
ing cross section. For a transition to a different umbrella 
state, or to a state with different j, we need potential 
energy coefficients with I > 1. Since viq and V30 are odd 
in p. they cause transitions from the 1 1 1 — ) state to the 
|21+) and the |11+) state for example. From Fig. HI we 
sec that the inelastic cross sections to these two states are 
indeed approximately equally large. The i^o term causes 
transitions to the 1 21 — ) state, and because this expan- 
sion coefficient is relatively large, the corresponding cross 
section is large as well. Finally, in order to change k in 
the collision, we need potential terms with to 7^ 0, of 
which the first two are the V33 and the U43 coefficients. 

The i>33 coefficient causes > — transitions and the 1143 

coefficient causes — — > + transitions. From Fig. [3] we 
see that the V33 coefficient is rather large, explaining the 
large cross sections to the |22— ) state, while the V43 coef- 
ficient is small, explaining the small cross sections to the 
|22+) state. 

In Fig. 1101 we show the integral cross sections for in- 
elastic scattering into the |32±) and the |31±) states, 
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FIG. 10: Inelastic integral cross sections for scattering into the 
|32±) and the |31±) states as a function of collision energy. 
At higher collision energies the |44±) channels for the para 
ammonia molecules open. As a result, Feshbach resonances 
are observed. 



for a collision energy ranging from the energies at which 
these channels open, up to 130 cm" 1 . These small 
cross sections will be harder to observe experimentally. 
However, if these channels can be detected, they allow 
for the observation of pronounced shape resonances at 
higher collision energies. The cross sections for scatter- 
ing into the |32±) states also give rise to Feshbach reso- 
nances with quasi-bound states of |31±) character. Even 
stronger Feshbach resonances at higher collision energies 
between about 120 and 125 cm -1 are found in the cross 
sections for scattering into the |31— ) state. These Fesh- 
bach resonances are caused by the opening of the |44±) 
channels at a collision energy of about 125 cm -1 . 

Finally, we also studied differential cross sections, 
where we looked in particular at the behavior of the dif- 
ferential cross sections as a function of energy close to 
resonance. In general, we found that the differential cross 
sections can change rapidly and dramatically close to res- 
onance. This behavior is illustrated by Fig. [TT] Here, we 
see in the upper two panels the differential cross sections 
for elastic scattering at collision energies of 0.87 and 1.86 
cm -1 . For the first of these energies there is no resonance, 
while for the second there is a shape resonance. For the 
off-resonance case we find that, apart from diffraction os- 
cillations, there is only a forward scattering peak. On res- 
onance there is also a strong backward peak. The lower 
four panels show the differential cross sections for inelas- 
tic scattering to the |11+) state at collision energies of 
20, 24.36, 30 and 37.28 cm -1 . At 20 and 30 cm" 1 , which 
are shown in the two lower plots on the left, there is no 
resonance and the differential cross sections look rather 
similar to the upper left one, giving predominantly rise 
to forward scattering. At 24.36 and 37.28 cm -1 , which 



are shown in the two lower plots on the right, there is 
a Feshbach resonance present, and as a result the dif- 
ferential cross sections look very differently, giving again 
rise to significant backscattering. In general, the precise 
structure of the differential cross section depends on var- 
ious aspects such as the lifetime and the rotational state 
of the intermediate collision complex. As a result, it is 
expected that the differential cross sections show clear 
changes near a resonance, but the precise way in which 
they change is hard to predict and can be very different 
for different resonances, as is also seen in Fig. [11] 



VI. DISCUSSION AND CONCLUSION 

In this article, we have theoretically studied cold colli- 
sions of NH3 molecules with He atoms, where we looked 
in detail at shape and Feshbach scattering resonances. 
Prior to collision, we considered the ammonia molecules 
to be in their antisymmetric umbrella state with angu- 
lar momentum j = 1 and projection k = 1, which is a 
suitable state for Stark deceleration. We calculated state- 
to-state integral and differential cross sections for colli- 
sion energies ranging from 10~ 4 cm -1 to 130 cm" 1 , using 
fully converged quantum close-coupling calculations. We 
treated the umbrella motion of the ammonia molecule 
by solving the corresponding Hamiltonian in curvilinear 
coordinates and taking the resulting first four vibration- 
tunneling states exactly into account. We call this the 
elaborate treatment. We also used a common model for 
the umbrella motion which approximates the umbrella 
tunneling states as an even and odd combination of the 
two possible rigid equilibrium structures for ammonia. 
This we call the model treatment. 

To describe the interaction between the NH3 molecules 
and the He atoms accurately, we presented a new four- 
dimensional potential energy surface, based on a high- 
quality fit of 4180 ab initio points. In the short range we 
used the explicitly correlated CCSD(T)-F12 method with 
an AVTZ basis set including midbond functions, and we 
showed that this approach leads to excellent results in 
the short range. In the long range we used the CCSD(T) 
method with an AVQZ basis but without F12, since we 
found that the electron density fitting that accompanies 
the F12 treatment does not exactly preserve charge neu- 
trality of the monomers and eventually leads to a dom- 
inant 1/R dependence of the potential at very large R 
values. Our potential has a well depth D e = 35.08 cm" 1 , 
which is to be compared with the well depth of 33.46 
cm" 1 for the potential of Hodges and Wheatley (29[. Al- 
though this difference is not very large, we found that 
small differences in the potential can have profound con- 
sequences for the observed resonance structures at low 
scattering energies. 

We studied all open collision channels for para ammo- 
nia up to j '■ — 3 and in all these channels we found pro- 
nounced shape resonances right after the opening of these 
channels, caused by quasi-bound states in the incoming 
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FIG. 11: Differential cross sections as a function of the zenith angle 9 for collisions between NH3 and He at various collision 
energies. Upper two panels: differential cross sections for elastic scattering at a collision energy of 0.87 cm" 1 (left) and 1.86 
cm" 1 (right). In the first case there is no resonance, the second case corresponds to a shape resonance. Middle two panels: 
differential cross sections for inelastic scattering into the |11+) state at a collision energy of 20 cm" 1 (left panel, no resonance) 
and 24.36 cm -1 (right panel, Feshbach resonance). Lower two panels: the same but for a collision energy of 30 cm -1 (left 
panel, no resonance) and 37.28 cm -1 (right panel, Feshbach resonance). 



and outgoing channels. We also found Feshbach reso- 
nances that are particularly strong for the outgoing |11+) 
channel at collision energies of about 25 cm -1 caused 
by intermediate |22±) states, and at collision energies of 
about 35 cm -1 caused by intermediate |21±) states. Due 
to the large cross section of these inelastic resonances, 
namely more than 1 A 2 , they seem to be a good candidate 
for experimental observation. Also in the |31— ) channel 
at collision energies of about 120 cm -1 relatively strong 
Feshbach resonances were seen that are due to interme- 
diate |44±) states. We analyzed the observed resonant 
structures in detail by looking at the corresponding scat- 
tering wavefunctions, phase shifts and lifetimes. We also 
investigated the validity of using the model treatment for 
the ammonia umbrella motion in describing low-energy 



scattering resonances. We found that the model per- 
forms qualitatively very well, but on a quantitative level 
some resonance peaks are somewhat overestimated com- 
pared to the elaborate treatment. However, considering 
the sensitivity of these resonances to the interaction po- 
tential, for which even state-of-the-art ab initio methods 
still lead to uncertainties on the order of a percent, the 
model treatment seems adequate in treating low-energy 
resonant scattering, especially in cases when the elabo- 
rate treatment becomes computationally too expensive. 

The calculated integral cross sections at low collision 
energies can be measured using Stark-decelerated molec- 
ular beams. The NH 3 molecule, and its isotopologuc 
ND3, arc amenable to the Stark deceleration technique, 
and have been employed frequently in deceleration ex- 
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FIG. 12: Convoluted inversion inelastic integral cross sections 
as a function of the mean collision energy for NH3-He scatter- 
ing. The initial state of the NH3 molecules is the 11 — ) state 
and the final state is the |11+) state. The figure is similar to 
Fig. only now we have assumed a Gaussian collision energy 
distribution for the colliding particles to simulate more realis- 
tically what would be observed with present day experimental 
technology. The blue curve corresponds to a full width at half 
maximum (FWHM) of 1 cm^ 1 , and the red curve of 3 cm^ 1 . 



pcriments |53[ . The Stark decelerator provides a beam of 
ammonia molecules, state-selected in the upper inversion 
component of the j = k = 1 level, with a velocity that is 
tunable between standstill and high velocities [54|. In a 
crossed beam experiment, the Stark decelerated ammo- 
nia molecules can be collided with an atomic beam of he- 
lium. In an optimized geometry, the two beams collide at 
a small beam intersection angle. An intersection angle of 
less than 90° reduces the attainable collision energy and 
improves the collision energy resolution of the experiment 
|55| . As shown in section |Vj there are a number of scat- 
tering channels with pronounced shape and/or Fcshbach 
resonances. The most promising prospects for the exper- 
imental observation of resonant features are found in the 
channels |11+) <- |11 — > (see Fig.©, |22-) <- 1 11 — ) , and 
|21-) <- |11-) (see Fig.©. 

The resonant structures are found at collision energies 
in the 20 - 50 cm -1 range. To simulate what would be 
observed in a molecular beam scattering experiment, we 
convoluted the integral scattering cross section for the 
1 11+) <— 1 11 — ) inversion inelastic channel with Gaussian 
collision energy distributions having both a 1 cm -1 and 
a 3 cm -1 full width at half maximum (FWHM). In the 
considered scattering channel bunches of Feshbach reso- 
nances are observed that are caused by the opening of 
the |22±) and |21±) channels, as seen in Fig. The re- 
sult of the convolutions are shown in Fig. From this 
figure we conclude that the details of the dense resonance 
structures in Fig. [3] can only be resolved when an exper- 



imental collision energy spread that is much less than 1 
enr 1 can be achieved. For an experimental resolution of 
3 cm -1 , however, the bunch of scattering resonances can 
still be discerned from the background inelastic signal by 
measuring the inelastic cross section as a function of col- 
lision energy. Figure [T^] shows that such a measurement 
would lead to a clear enhancement of the inelastic signal 
by more than a factor of two at the energies where the 
bunch of resonances is located. 

To estimate the feasibility of obtaining collision en- 
ergy resolutions in this range with current experimental 
technology, we assume an experiment in which a Stark- 
decelerated packet of NH3 molecules collides with a con- 
ventional beam of He atoms at a beam intersection angle 
of 45°. We assume that the He atom beam is produced 
using a cryogenic source that is maintained at a temper- 
ature of about 50 Kelvin, resulting in a He atom velocity 
of 550 m/s. The relevant range of collision energies is 
then obtained when the velocity of the NH3 molecules 
is tuned between 370 and 850 m/s. This is well within 
the range of state-of-the-art Stark deceleration molec- 
ular beam machines [Htjj . We further assume velocity 
spreads of 10 m/s and 55 m/s for the NH3 molecules and 
He atoms, respectively, and a spread in the beam inter- 
section angle of 40 mrad due to the divergence of both 
beams. These are values that can realistically be ob- 
tained with current experimental techniques. With these 
parameters, we expect an optimum in the collision en- 
ergy resolution to occur at a collision energy of 30 cm -1 , 
i.e., at the center of the relevant collision energy range. 
This maximum accuracy amounts to a spread of 3.1 cm -1 
(FWHM), while the collision energy spread increases to 
approximately 4 cm -1 both for collision energies down to 
20 cm -1 and collision energies up to 50 cm -1 . These en- 
ergy resolutions will not yet allow for the observation of 
single scattering resonances in the NH3-He system, but 
they will certainly enable to observe the enhancement of 
about a factor of two in the inelastic cross section as a 
function of collision energy shown in Fig. 1121 revealing 
the combined effect of the underlying bunch of Feshbach 
resonances. 

An alternative and complementary approach to study 
scattering resonances is to measure differential cross sec- 
tions. Referring back to Fig. [TTJ dramatic changes in 
the differential cross section can occur at collision ener- 
gies where a resonance is observed. Feshbach resonances 
that give rise to strong backward scattering can be de- 
tected by measuring selectively the scattered flux in the 
backward direction. A similar approach has been used re- 
cently to measure partial- wave resolved resonances in the 
collision energy dependence of reactions between F atoms 
and HD molecules @. For inelastic scattering between 
NH3 molecules and He atoms, differential cross sections 
are measured most conveniently using the velocity map 
imaging (VMI) technique [57| . To experimentally resolve 
the angular dependence of the differential cross sections, 
large recoil velocities of the scattered molecules are ad- 
vantageous. For the |11+) <— |H— ) channel, the recoil 
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velocity of the scattered NH3 molecules in the center of 
mass frame amounts to about 100 m/s at the most rel- 
evant collision energies. This is well within the range of 
velocities that can be imaged using current VMI tech- 
niques, offering interesting prospects to study the behav- 
ior of molecular scattering resonances. 

At the Fritz-Habcr-Institute in Berlin, Germany, and 
the Radboud University Nijmegen, The Netherlands, we 
have embarked on an experimental program to study 
scattering resonances in both the integral and differential 
cross sections using Stark-decelerated molecular beams. 
It is the hope that the experimental study of these reso- 
nances will test our theoretical understanding of molec- 
ular interactions with unprecedented accuracy, and also 
will contribute to an enhancement of our ability to con- 
trol the way in which molecules collide. 



Appendix A: Coordinates and symmetry 

In this Appendix we consider in more detail the vari- 
ous coordinates used to describe the NH3-He system and 
how these coordinates transform under various symmetry 
operations that commute with the Hamiltonian. These 
symmetry operations form a group generated by the per- 
mutations of the hydrogen atoms in NH3 and the op- 
erator for inversion in space E* . The location of the 
center-of-mass of the NH3-He dimer is given by the vec- 
tor Q, while the locations of the nuclei with respect to 
Q are given by the vectors Hi, H2, H3 for the three 
H atoms, by the vector N for the N atom, and by He 
for the He atom. The center-of-mass of the ammonia 
molecule is given by the vector X. We start by introduc- 
ing a orthonormal, right-handed space-fixed ('sf') frame 
centered at the center-of-mass of the dimer Q. We also 
make the convention that a superscript denotes the frame 
in which the coordinates of a vector are expressed. So 
in the space-fixed frame an arbitrary vector P has the 
space-fixed coordinates P sf . Moreover, when no frame 
superscript is given, we do not specify the frame in which 
the coordinates of the vector P are evaluated. The space- 
fixed frame consists of three unit vectors, which are the 
columns of the matrix s = (s^, s y , s z ). In our notation s sf 
is the unit matrix. Moreover, in the space-fixed frame the 
inversion operator E* maps the position of any nucleus 
P sf onto the position reflected in the origin — P sf . 

Another useful frame, the dimer frame or body-fixed 
frame d = (A x ,A. yi A z ), is obtained by performing two 
rotations to the space-fixed frame in order to align the 
d z axis of the dimer frame along the vector R = He — 
X, which points from the center-of-mass of the ammonia 
molecule to the helium atom. We have in space-fixed 
coordinates that 



' i?cos0 sf sin6» sf 
i?sin0 sf sin0 sf 
^ Rcos9 sl 



(Al) 



so that (j) si and 9 st are the azimuth and zenith angle of the 
vector R in the space-fixed frame. The body-fixed frame 
is then defined in the following way. For any vector P, we 
have that P = s • P sf = d • P bf with P sf = TV\\ • P bf and 



are 



7t« = 7l z ((t) Bt )7l y (e Bt ), where TV v {6 st ) and 1Z Z 
the usual rotation matrices for rotation about the y axis 
and the z axis, respectively, with the convention that 
TV z ((t> si ) 12 = -sin</) sf and 7^(6> sf ) 3 i = — sin 6» sf . Note 
that as a result we find for the coordinates of the vector R 
in the body-fixed frame that (TV 



bi) 



R bf resulting 

in R x l = Ry l = and i? bf = i?, as was required. When 

we apply the E* operator to the complex, the coordinates 
R sf get inverted, so that E*: R sf — > — R sf . As a result, 
the angles in Eq. (|A1[) are changed according to </> sf — >• 
4> sl + tt and (9 sf — > tt — 9 sl . Note that the dimer frame is 
invariant under permutations of the hydrogen atoms. 

The third useful frame is called the monomer frame 
and it is located at the center-of-mass of the ammonia 
molecule X. The frame is spanned by the vectors 



w x = 2Hi H 2 H 3 , 

Vy = H 2 — H 3 , 
V Z = V X X Vy. 



(A2) 



Since the ammonia molecule keeps it threefold symme- 
try, the above frame is orthogonal. The monomer frame 
can consequently be made orthonormal by dividing the 
vectors in Eq. (|A2[) by their length, resulting in the 
right-handed monomer frame denoted by f . The rota- 
tion matrix that expresses the monomer frame axes in 
space-fixed coordinates is given in terms of the three 
Euler angles Q si — (a si , /3 st , 7 sf ), resulting in TV s mt = 
TV z (a sf )TVy(p si )TV z (Y { )- When the E* operator is ap- 
plied to the complex, then both the x axis and the y axis 
of the monomer frame are reversed in the space-fixed 
frame, which leaves the z axis in place. As a result, the 
E* operator has the following effect on the Euler angles, 
E*: a sf -> a sf , /3 sf -> /3 sf , and 7 sf ->■ 7 sf + tt. Another 
angle that is important in our treatment of the ammonia- 
helium complex is the inversion or umbrella angle p, de- 
fined as the angle between the z axis of the monomer 
frame and the vector pointing from the N atom to one 
of the H atoms. So, for p = tt/2, ammonia has a planar 
geometry. We just showed that E* leaves the z axis in 
place, while the coordinates of the nuclei change sign. As 
a result, we have that E*: p — > n — p. 

When we permute the hydrogen nuclei, the space-fixed 
coordinates of the monomer frame axes are interchanged. 
For example, when we interchange H 2 and H3 we find 
that the monomer y axis and also the z axis are re- 
versed in the space-fixed frame. As a result, the corre- 
sponding Euler angles transform as (23): a sf — > tt + a sl , 
/3 sl -> tt - /3 sf , and 7 sf -> - 7 sf . Since the (23) op- 
eration inverts the monomer z axis, we also find that 
(23): p — > 7r — p. In Table [TTI we show the transforma- 
tion properties of the various angles that describe the 
ammonia-helium complex when symmetry operations of 
the permution- inversion group D^h (M) are applied to the 
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TABLE II: Transformation properties of the coordinates under symmetry operations. 



Operation Effect on angular coordinates 



E 




e si 








p 




a sf 




7 sf 


a bf 




/3 bf 




7 bf 


(123) 












p 




a sf 




7 sf -2tt/3 


a bf 




0bt 




-2tt/3 


(23)* 


7T 


- e sl 


4> st + 


it 




p 




a sf + 7r 


7T-/3 sf 


7T -7 sf 


7T-a bf 




P hi 




-7 M 


E* 


7T 


-6 sC 


<t> st + 


it 






p 


a 


I3 st 


7 Sf +7T 


bf 

—a 


TT 


-/3 bf 




7 bf 


(123)* 


7T 


_0Sf 


<f> si + 








p 


a 




7 sf + t/3 


bf 

—a 


71" 


_ /3 m 




-2tt/3 


(23) 




9 si 






IX 




p 


a Bt + % 


7T-/3 sf 


-7 sf 


a bf + 7T 


TT 


-/3 bf 




-7 W 



TABLE III: Effect of the symmetry operations on the angular basis functions. 



Operation Body-fixed Space-fixed 

E \jkKJMv ± ) IjkLJMv^ 

(123) e 2nlk/3 \jkKJMv ± ) e 27Tlk/3 \jkLJMv ± ) 

(23)* ±(-l) J+k \j-k-KJMv ± ) ±(-l) j+k+L \j-kLJMv ± ) 

E* ±{-l) J+J+k \jk-KJMv ± ) ±{~l) L+k \jkLJMv ± ) 

(123)* ±(-l) J+3+k e 27Ilk/3 \jk-KJMv ± ) ±e 2nik/3 \jkLJMv ± ) 

(23) (-iy\j-kKJMv ± ) (-l) J \j-kLJMv ± ) 



complex. These transformation properties are useful in 
determining the transformation properties of the angular 
basis functions. 

The monomer frame can also be obtained by a rota- 
tion from the body-fixed dimer frame rather than the 
space-fixed frame, defining the body-fixed Euler angles 
C bf = (a bf , /3 bf , 7 bf ). The rotation matrix that expresses 
the monomer frame axes in body-fixed coordinates is 
given by 7?. bf f = n z (a hi )n y (P hf )n z {-/ hf ). When the 
hydrogen atoms are permuted, the body-fixed frame is 
unchanged, and as a result the body-fixed Euler an- 
gles transform in precisely the same way as the space- 
fixed ones. However, when E* is applied, not only the 
monomer frame axes in space-fixed coordinates change, 
but also the dimer frame axes. As a result, the body- 
fixed Euler angles transform somewhat differently than 
the space-fixed angles, as seen in Table HH 

Having determined the effect of the various symmetry 
operations on the coordinates that describe the NHrHe 
complex, we can also find out the corresponding effect on 
the angular basis functions of Eqs. (J5J) and ([7]) by using 
the transformation properties of the Wigner d-functions. 
Moreover, we have that E*(fy^(p) = ±<j)f(p). The effect 
of the various symmetry operations on the angular basis 
functions is shown in Tabic [TTTT 

As a result, we are now able to construct the symme- 
try adapted bases sets for both the body-fixed and the 
space-fixed case. To this end, it is most convenient to 
start by discussing the symmetry group C% v (M) with ir- 
reps A\, A2 and E. The Cj, v {M) group is generated by 
the operations E, (123) and (23)*. Using this group im- 
plies that we consider the ammonia molecule as a rigid 
rotor without umbrella motion. Then, for the body-fixed 
case with k = K = 0, we conclude from Table Mil that 
the state \j00JM) is of A\ symmetry when J is even, 
while it is of A2 symmetry when J is odd. When ei- 



ther k or K is nonzero, we have for k = (mod 3) that 
(\jkKJM) + {-l) J+k \j-k-KJM))/2^ 2 is of A 1 sym- 
metry, while (\jkKJM) - (-l) J+k \j-k-KJM))/2 1 / 2 is 
of A2 symmetry. Finally, we have for k 7^ (mod 3), 
that two-dimensional E irreps are spanned by the states 
(\jkKJM),\j-k-KJM)). 

For the space-fixed case and considering rigid ammo- 
nia, we conclude from Table HTT1 that the state \jOLJM) 
is of A\ symmetry when j + L is even, while it is 
of Ai symmetry when j + L is odd. When k is 
nonzero and a multiple of 3, we have that (\jkLJM) + 
(-iy +L+k \j-kLJM))/2 1 / 2 is of Ax symmetry, while 
(\jkLJM)-(-iy +L+k \j-kLJM))/2 1 / 2 is of A 2 symme- 
try. For k 7^ (mod 3), we have that two-dimensional E 
irreps are spanned by the states {\jkLJM), \j—kLJM)}. 

Finally, the complete symmetry adapted basis is ob- 
tained by considering the full D^h (M) symmetry group of 
the nonrigid ammonia molecule. The functions adapted 
to the irreps of D^h (M) are obtained from those adapted 
to the C^ v (M) irreps by using 

(E + E*)\Ai) = \A[), (E - E*)\Ai) = \A%), 
(E + E*)\A 2 ) = \A' 2 ), (E - E*)\A 2 ) = \A'{), 
(E + E*)\E) = \E'), (E-E*)\E) = \E"). (A3) 

Further we note that the v + umbrella functions belong to 
the (E + E*) projection, while the v~ umbrella functions 
belong to the {E — E*) projection. 
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